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Abstract 

Using a sample of 2.59 x 10 7 ip(2S) decays collected by the CLEO-c detector, we present results 
of a study of Xcj (.7=0,1,2) decays into baryon-antibaryon final states. We present the world's 
most precise measurements of the Xcj — » PP and Xcj —* AA branching fractions, and the first 
measurements of XcO decays to other hyperons. These results illuminate the decay mechanism of 
the Xc states. 

PACS numbers: 13.25.Gv, 14.40. Gx 
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In the standard quark model, the Xcj {J = 0,1,2), mesons are cc states in an L = 1 
configuration. The Xcj mesons are not produced directly in e + e~ annihilations. However, 
the large branching fractions of i[>(2S) — > XcJl make e + e~ collisions at the ip(2S) energy a 
very clean environment for Xcj investigation. 

The available data on the decays of the Xcj mesons into baryon-antibaryon pairs has so 
far been very limited. The easiest of these final states to detect and measure is pp [l[ ■ The 
partial width for Xco ~^Vp was originally predicted to be zero in some models due to the 
Helicity Selection Rule |2j. However, this rule has long been known to be strongly violated. 
More recent work has concentrated on the importance of the Color Octet Mechanism (COM), 
which treats the Xc states as more than just pure qq states and incorporates octet operators in 
the transition matrix elements to a given final state in order to calculate two-body exclusive 
decay rates [3J. In particular, Wong [4j used the COM to explain the high rate of Xcj PP 
and made predictions for Xcj —> AA. However, these predictions fell well below the low- 
statistics measurements from BES [B| that imply B(xc ~^ AA)/i5(x c ~^ PP) ~ 2 to 4 for all 
three Xc states. It has since been postulated that such large ratios can be explained without 
using the COM, and, instead including a more detailed quark model of the daughter products 
0. However, the resulting predictions depend greatly on the details of this model, and it is 
clear that more experimental input is needed. In this paper, we analyze a large sample of 
i/j(2S) decays and present results on two-body decays of the Xcj mesons into pp, AA, S°S , 



The data were taken by the CLEO-c detector [7| operating at the Cornell Electron Storage 
Ring with e + e~ collisions at a center of mass energy corresponding to the ip{2S) mass of 
3.686 GeV/c 2 . The data correspond to an integrated luminosity of 56.3 pb -1 , and the total 
number of i[>(2S) events is calculated as 2.59 x 10 7 , determined according to the method 
described in [8J. 

Photons were detected using the Csl crystal calorimeter which has an energy resolu- 
tion of 2.2% at 1 GeV, and 5% at 100 MeV. To discriminate protons from kaons and pions, 
we combined specific ionizations (dE/dx) measured in the drift chamber and log-likelihoods 
obtained from the ring-imaging Cerenkov detector (RICH) [Io| to form a log-likelihood dif- 
ference: £(p - 7r) = Crich(p) - Crich{k) + aj E/dx (p) - o\ Eldx {n), where negative C(p - vr) 
implies the particle is more likely to be proton than a pion. For all protons in the events we 
require C(jp — ir) < and C(p — K) < 0. This is a very efficient requirement. 

We reconstruct the hyperons in the following decay modes: A — > pn~ (branching fraction 



63.9%) Q , E+ -►jwr°(51.6%),E° -> A 7 (100%),~- -> Att-(99.9%_K and S° -> Att°(99.5%). 



Our hyperon detection follows the technique explained elsewhere [111]. Briefly, to reconstruct 
A candidates, proton candidates are combined with charged tracks that are assumed to be 
pions. The pir combination is required to be within 10 MeV of the known A mass and then 
is kinematically constrained to that value. Similarly, S~ candidates are built from these A 
candidates with the addition of another appropriately charged track assumed to be a pion. 
The An vertex was required to be closer to the beamspot than the A decay point. The 
S° candidates were formed from the combination of A candidates and a cluster of greater 
than 50 MeV energy detected in the crystal calorimeter, not matched to the trajectory of a 
charged track, and consistent in shape with that expected from a photon. The S + and S° 
reconstruction is complicated by the fact that we cannot use the beamspot for the point of 
origin of the photons. A kinematic fit is made to the hypothesis that the parent hyperon 
originated at the beamspot, and decayed after a positive path- length at a point taken to 
be the origin of the it — > 77 decay. A requirement was placed on the x 2 °f the fit to 



E+E+, 



and H°S . 
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this topology, which includes the fit to the ir° mass from the newly found decay vertex. In 
all cases, hyperon candidates within 3a of their nominal masses are considered for further 
analysis, and their four-momenta are then constrained to the nominal hyperon mass. These 
kinematic constraints were sufficient to ensure that cross-feed background from real Xcj 
decays, for instance AAtt + ti^ in the sample, was negligible. 

For events with two distinct baryon candidates, we combine the candidates into a Xc 
candidate. At this stage of the analysis, the invariant mass resolution of the Xc is around 
15 MeV/c 2 . We then search for any unused photon in the event and add that to the Xc 
candidate to form a ip{2S) candidate. This ip(2S) is then kinematically constrained to the 
four- momentum of the beam, the energy of which is calculated using the known ip(2S) 
mass. The momentum is non-zero due to the finite crossing angle (~ 3 mrad per beam) 
in CESR. To make our final selection, we require the ip{2S) candidate to have a x 2 °f l ess 
than 25 for the four degrees of freedom for this fit; this requirement rejects most background 
combinations. This kinematic fit greatly improves the mass resolution of the Xc candidate. 

To study the efficiency and resolutions, we generated Monte Carlo samples for each Xc 
into each final state using a GEANT-based detector simulation [12]. The simulated events 
have an angular distribution of (1 + a cos 2 9), where 9 is the radiated photon angle relative to 
the positron beam direction, and a = 1, -1/3, and 1/13 for the Xco, Xa, and Xc2 respectively, 
in accordance with expectations for an El transition. The mass resolution and efficiencies 
are shown in Table I. The resolutions are approximated by single Gaussian signal functions. 
The efficiencies shown include all the relevant branching fractions. [l[ 

TABLE I: Efficiencies (in %) obtained from analysis of Monte Carlo generated events, and yields 
found in the data sample. 



Mode 


XcO 


Xcl 


Xc2 




Yield 


Efficiency (%) 


Yield 


Efficiency (%) 


Yield 


Efficiency(%) 


PP 


383 ± 22 


62.4 


141 ± 13 


66.6 


121 ± 12 


65.5 


AA 


131 ± 12 


16.2 


46.0 ± 7.2 


17.1 


71.0 ±9.2 


17.3 


S°E° 


78 ± 10 


4.1 


3.8 ± 2.5 


4.0 


7.5 ±3.4 


4.0 


£+£+ 


39 ±7 


5.2 


4.3 ± 2.3 


5.0 


4.0 ±2.3 


4.7 




95 ±11 


7.7 


16.4 ±4.3 


8.2 


29 ±5 


8.4 




23.3 ±4.9 


2.9 


1.7 ± 1.4 


2.9 


2.9 ± 1.7 


2.9 



The final invariant mass distributions are shown in Fig. 1. These plots are each fit with 
three signal shapes comprising Breit-Wigner functions convolved with Gaussian resolutions, 
together with a constant background term. The masses and widths of the Breit-Wigner 
functions were fixed according to the current averages [l[, and the widths of the Gaussian 
resolution functions were fixed at the values found from Monte Carlo simulation (ranging 
from 3.6-5.1 MeV/c 2 depending on the spin of the Xc and the decay mode). The yields from 
these fits are tabulated in Table I. 

To convert the yields to branching fractions, we divide by the product of the number 
of ip(2S) events in the data sample, the detector efficiency, and the branching fractions for 
i/j(2S) into Xcj- For the last factor we use the CLEO measurements of B(if)(2S) — > 7Xco) = 
9.22 ± 0.11 ± 0.46%, B(ip(2S) -> 7Xd) = 9.07 ± 0.11 ± 0.54%, and B{ip{2S) -> 7x c2 ) = 
9.33 ± 0.14 ± 0.61% [8J. The results are tabulated in Table II. 

We consider systematic uncertainties from many different sources. All modes have a 2% 



4 



uncertainty from the total number of ij)(2S) decays [8[]. The requirement on the \ 2 of the 
constraint to the beam four-momentum has been checked by changing the cut and noting the 
change in the yield in these, and other similar decay modes. Based on this study we place 
a systematic uncertainty of 2.5% on the efficiency of this requirement. The uncertainties 
due to track reconstruction are 0.3% per charged track. The limited Monte Carlo statistics 
introduces an uncertainty that is always a small fraction of the statistical uncertainty in the 
data. Using comparison of data and Monte Carlo simulation of hyperon and anti-hyperon 
yields from the ip(2S), we checked our modeling of the hyperon selection efficiency. The 
assigned systematic uncertainty arising from this study was up to 3% per hyperon. The 
systematic uncertainty due to the photon detection and shower-shape criteria is set at 2% 
per photon. In the case of the Xci decaying into two spin one-half particles, the two daughters 
can have their spins either parallel or antiparallel, and in the Xc2 case there are even more 
possibilities of combinations of intrinsic spins and relative angular momentum. These helicity 
correlations are not well known in the case of decays into baryons, and this introduces a 
small uncertainty in the modeling of the efficiencies. We investigated the effects of helicity 
amplitudes on our efficiency by generating Monte Carlo with a variety of different helicities 
and found small variations. From this study, we assign a 1% uncertainty in the efficiency 
of the Xci and 2.5% of the Xc2- When calculating the final branching fractions, we add the 
above systematic uncertainties in quadrature. The uncertainty due to the ij)(2S) — > jXc 
branching fractions is kept separate and quoted as a second systematic uncertainty. 

For evaluating the limits in the cases where there is no significant signal, we take the 
probability density function and convolve this with Gaussian systematic uncertainties. We 
then find the branching fraction that includes 90% of the total area. 

TABLE II: Branching fraction results (in units of 10~ 5 ) for each decay mode. The uncertainties 
are statistical, systematic due to this measurement, and systematic due to the ip{2S) — > XcJl rate, 
respectively. The limits on the branching fractions include all systematic uncertainties. 



Mode 




XcO 


Xcl 


Xc2 


PP 


This Work 


25.7 ±1.5 ±1.5 ±1.3 


9.0 ±0.8 ±0.4 ±0.5 


7.7 ±0.8 ±0.4 ±0.5 




PDG 


22.5 ±2.7 


7.2 ± 1.3 


6.8 ±0.7 


AA 


This Work 


33.8 ± 3.6 ± 2.2 ± 1.7 11.6 ± 1.8 ± 0.7 ± 0.7 17.0 ± 2.2 ± 1.1 ± 1.1 




PDG 


47.0 ± 16.0 


26.0 ± 12.0 


34.0 ± 17.0 




This Work 


44.1 ±5.6 ±4.2 ±2.2 


< 4.4 


< 7.5 




PDG 








£+£+ 


This Work 


32.5 ±5.7 ±4.0 ±1.7 


< 6.5 


< 6.7 




PDG 










This Work 


51.4 ±6.0 ±3.9 ±2.6 


8.6 ±2.2 ±0.6 ±0.5 


14.5 ±3.0 ±1.2 ±0.9 




PDG 


< 103 a 


< 34 


< 37 




This Work 


33.4 ±7.0 ±4.5 ±1.7 


< 6.0 


< 10.6 




PDG 









a The BES central value [13j for this measurement is (53 ± 27 ± 9) x 10 5 , in good agreement with this 
work. 

In summary, we measure branching fractions for XcO decays into pp, AA, 

and S + S+. For Xci an d Xc2 we find significant signals and measure branching frac- 
tions into the first three of the above decay modes. Upper limits on branching fractions 
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are obtained for the remainder of the modes. In the case of Xcj —>■ PP and Xcj —* AA, 
these measurements are the most precise to date; in the other modes they represent first 
measurements. Our values of the branching fractions for AA are well below those reported 
by BES, but confirm the trend that the branching fractions into A A are higher than those 
for pp. The fact that the Xco branching fractions into EE and 55 are all greater than that 
°f Xco ~ > PP, a trend not mirrored in the Xci and Xc2 decays, is not in agreement with naive 
expectations for the decay of an SU(3) singlet. 
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luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
partment of Energy, the Natural Sciences and Engineering Research Council of Canada, and 
the U.K. Science and Technology Facilities Council. 
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